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ABSTRACT 



Context. Methanol masers at 6.7 GHz have been found exclusively towards high-mass star forming regions. Recently, some Class 
protostars have been found to display conditions similar to what are found in hot cores that are associated with massive star formation. 
These hot corino sources have densities, gas temperatures, and methanol abundances that are adequate for exciting strong 6.7 GHz 
maser emission. 

Aims. This raises the question of whether 6.7 GHz methanol masers can be found in both hot corinos and massive star forming 
regions, and if not, whether thermal methanol emission can be detected. 

Methods. We searched for the 6.7 GHz methanol line towards five hot corino sources in the Perseus region using the Arecibo radio 
telescope. To constrain the excitation conditions of methanol, we observed thermal submillimeter lines of methanol in the NGC1333- 
IRAS 4 region with the APEX telescope. 

Results. We did not detect 6.7 GHz emission in any of the sources, but found absorption against the cosmic microwave background in 
NGC1333-IRAS 4A and NGC1333-IRAS 4B. Using a large velocity gradient analysis, we modeled the excitation of methanol over a 
wide range of physical parameters, and verify that the 6.7 GHz line is indeed strongly anti-inverted for densities lower than 10^ cm"''. 
We used the submillimeter observations of methanol to verify the predictions of our model for IRAS 4A by comparison with other 
CH3OH transitions. Our results indicate that the methanol observations from the APEX and Arecibo telescopes are consistent with 
dense (n ~ 10* cm^^), cold {T ~ 15 - 30 K) gas. 

Conclusions. The lack of maser emission in hot corinos and low-mass protostellar objects in general may be due to densities that are 
much higher than the quenching density in the region where the radiation field is conducive to maser pumping. 

Key words. Masers - Stars: low-mass 



' 1. Introduction 

■ The 5i - 60 A^ line of methanol at 6.7 GHz is the strongest 
[ of Class II methanol masers. 6.7 GHz masers are unique 
. compared to OH and H2O masers in that they are associ- 
' ated exclusively with high-mass star forming regions. For ex- 
ample, targeted searches towards low-mass young stellar ob- 
jects ha ve not resulted in any detections of 6.7 GHz methano l 
masers (! Minieretal.ll200l iBourke et aLllIOOl IXu etal.ll2008h . 
Theoretical models suggest that the maser emission switches on 
at gas densities greater than ~ 10"* cm^^and methanol fractional 
abundance greater than ~ 10"^, along with the requirement of a 
far-infrared radiation field to excite metha nol to torsionally ex- 
cited sta tes as par t of th e pumping cycle dSobolev et alJll997t 
ICragg et al. 2002, l 2005h . These conditions are found in high- 
mass star forming regions where the far-infrared radiation field 
is generated by warm dust (To ~ 100 - 150 K), and high 
methanol abundances are achieved by evaporating the mantles of 
dust grains, either by the passage of a shock or through the radi- 
ation field of the massive young stellar object. The molecules re- 
leased from the grains further trigger f ormation of more com plex 
molecules in warm gas chemistry (e.g. lCharnlev et anil992h . 

Recently, there has been evidence of hot core chemistry in 
low-mass protostars. For example, complex molecules such as 
HCOOCH3, HCOOH, and CH3CN have been found in IRAS 



16293-24 22 ("Cazaux elalj l2003h and NGC1333-IRAS 4A 
jBottinelli et al. ,2004J along with high abun da nce of H2O, 
CHiO H, and H2CO (ICeccarem et al J 1200035 ISchoier et all 
'20021). The abundance of methanol in the compact, warm regions 
in the immediate proximity of these emb edded sources, kno wn 
as hot corinos, can be as high as 7 x 10~^ dMaret et alj|2005l and 
references therein). Moreover, the regions have warm tempera- 
ture (> 100 K) and high density (> 10'' cm~^; ICeccarelh et all 
12000 a'). This raises the question of whether hot corinos can 
excite maser emission in Class II methanol lines. A discov- 
ery of Class II methanol masers in hot corinos would force re- 
interpretation of the methanol maser data in the literature. In this 
paper, we report on the results of a search for 6.7 GHz methanol 
masers towards five hot corino sources in the Perseus region. 

2. Observations 

2.1. Arecibo Observations 

The search for 6.7 GHz maser emission towards hot corinos was 
caiTied out using the 305 m Arecibo radio telescop43 in October 
2006. The C-High receiver was used along with the "interim" 



' The Arecibo Observatory is part of the National Astronomy and 
Ionosphere Center, which is operated by Cornell University under a co- 
operative agreement with the National Science Foundation. 
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Table 1. Hot corinos observed at 6.7 GHz. The columns show 
source name, equatorial J2000 coordinates and the rm.s. noise 
in spectra. 



Source 


a 


S 


5 (rms) 




(J2000) 


(J2000) 


mjy 


L1448-N 


03 25 36.6 


+30 45 15 


1.2 


L1448-MM 


03 25 38.8 


+30 44 05 


1.3 


NGC1333-IRAS2 


03 28 55.4 


+31 14 35 


1.2 


NGC1333-IRAS 4A... 


03 29 10.3 


+31 13 31 


1.0 


NGC1333-IRAS 4B... 


03 29 12.0 


+31 13 09 


1.0 



correlator backend. Two units of the correlator were used with 
1.5625 MHz bandwidth to measure two orthogonal linear po- 
larizations with 9-level sampling and 1024 channels per polar- 
ization giving a velocity coverage of 140 km s"' with 0.28 km 
s ' resolution after Hanning smoothing. The adopted rest fre- 
quency for the transition was 6668.518 MHz; a more accurate 
measurement of 6668.519 MHz by Breckenridge & Kukolich 
gives an error of 0.04 km s ' in the velocity scale, which 
is well below the resolution of the spectra. The sources ob- 
served and their coordinates are indicated in Table 1. Sources 
NGC1333-IRAS 4A and NGC1333-IRAS 4B had on-source in- 
tegration times of 45 minutes each, while the remaining sources 
had on-source integration times of 30 minutes each. The data 
were taken in position switched mode. The "off" spectra were 
taken over the same elevation and azimuth range as for the "on" 
spectra to ensure cancellation of standing waves. During each 
observation, the pointing was checked using a continuum source 
and was found to be accurate within 10". The system tempera- 
ture was between 23 and 34 K, with approximately 20% differ- 
ence between the two polarizations. The half-power beamwidth 
(HPBW) of the telescope at this frequency is ~ 40". 

The data were reduced in IDL using procedures maintained 
by the observatory. The antenna temperature was derived using 
noise diodes, which was converted to flux density using an ele- 
vation and azimuth dependent gain curve for this frequency, the 
typical conversion factor being 5 K/Jy. The resulting rm.s. noise 
in the spectra after Hanning smoothing was 1 .0 mJy for IRAS 4A 
and IRAS 4B, and 1.2-1.3 mJy for the remaining sources. We 
did not carry out any baselining of the spectra. 

2.2. APEX Observations 

It is possible to carry out an excitation analysis to reproduce the 
observed 6.7 GHz feature using milli meter and submillimeter 
lines of methanol jLeurini et al.ll2004 ). In order to do this, we 
observed the NGC1333-IRAS 4 region in the (6jf - 5k) CH3OH 
band at ~ 290 GHz using the APEX 12 m telescopeQ. The obser- 
vations were carried out in November 2006 using the APEX2A 
double side band receiver dRisacher et al.l |2006|) and th e Fast 
Fourier Transform Spectrometer (FFTS; l KIeinet al.ll2006 ). The 
receiver was tuned to a frequency of 290.300 GHz in the lower 
sideband. One unit of the FFTS was centered at the tuning fre- 
quency, while the second unit was centered at 291.5 GHz. Both 
units had a bandwidth of 1 GHz with 16 384 channels, the ef- 
fective velocity resolution being 0.13 km s~' (since the noise in 
successive channels is correlated: iKlein et al.l l2b06). The point- 
ing was checked approximately every 1.5 hours in CO (J - 3-2) 

^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



on IK-Tau and was found to be accurate to within 2". The focus 
was verified using Saturn. The system temperatures were typ- 
ically 120-150 K in units of T^. The half -power beamwidth of 
the telescope at this frequency is around 21". A 1.5' x2. 5' region 
was mapped using the on-the-fly scanning mode with a step size 
of 6". 

The data were reduced using the CLASS software. A first 
order polynomial was fitted to the line-free channels and sub- 
tracted off the baseline. The data were converted to main beam 
temperature units (Tmb) ass uming a be am efficiency of 0.73 and 
a forward efficiency of 0.97 (iGtisten et al. 2006). In dividual lines 
were then identified using the XCLASS software (IComito et al.l 
l2005h . 

3. Results 

We did not detect 6.7 GHz maser emission in any of the sources, 
but detected prominent absorption in NGC1333-IRAS 4A 
and NGC1333-IRAS 4B, as shown in Fig. [T] IRAS 4A and 
IRAS 4B are separated by 31", which is less than the half- 
power beamwidth of 40". Thus, the measured spectrum towards 
IRAS 4A has some contribution from IRAS 4B and vice versa. 
Assuming a symmetric beam, if a is the contribution of the sec- 
ond source at the position of the first, then the measured spectra 
can be written as 

'^IRAS4A ~ IRAS4A ^ "'^ IRAS4B y'^' 
'^IRAS4B - '^i'rAS4B + '^'^IRAS4A (2) 

where S" and S'" indicate the actual and measured spectra re- 
spectively. Since the beam is a Gaussian to good approxima- 
tion at 31" from center, the parameter a is calculated as 0.19. 
Equations (1) and (2) can then be inverted to calculate the ac- 
tual spectra, the results being shown in Fig.|2] The contribution 
of IRAS 4B to the observed IRAS 4A spectrum is small (less 
than ~ 10%), while the contribution of IRAS 4A to the observed 
IRAS 4B spectrum is much more significant (~ 30%). 

For modeling the observed 6.7 GHz absorption, we 
smoothed the APEX data to the spatial scale of the Arecibo 
beam. The resulting spectrum at the IRAS 4A position is shown 
in Fig. [3] We were unable to obtain spectra from the smoothed 
data cube at the IRAS 4B position since the footprint of the 
Arecibo beam at this location extends beyond the eastern edge of 
the APEX map. Since the convolution will give rise to the same 
degree of blending between the two sources as in the Arecibo 
observation, one would ideally have to deconvolve the spectra 
of the two sources extracted from the smoothed data cube, using 
the same parameters as for the Arecibo data. In the absence of 
spectra from both IRAS 4A and IRAS 4B at the same spatial res- 
olution, it is not possible to deconvolve the two sources as was 
carried out for the Arecibo data. To quantify this contamination 
in the APEX data, we extracted a spectrum from IRAS 4A tak- 
ing care to exclude the IRAS 4B region during the spatial convo- 
lution, and compared it with the spectrum at the same position 
from the original smoothed datacube (which includes the con- 
tribution from IRAS 4B). Fig |4] shows the difference between 
the two spectra overlaid on the IRAS 4A spectrum from the 
smoothed datacube. From this analysis, we estimate the contri- 
bution of IRAS 4B to IRAS 4A in the APEX data to be at most 
3.5%. Henceforth, we will ignore the contamination of IRAS 4B 
on IRAS 4A. Since the contribution of IRAS 4A to the IRAS 4B 
spectrum is significant, we do not model IRAS 4B. 

Fig. |5] shows the 6.7 GHz absorption spectrum measured 
towards IRAS 4A and the 290.1 GHz thermal methanol emis- 
sion over the same velocity range. The similarity of the line 
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Fig. 1. The observed 6.7 GHz spectra showing absorption in Fig. 3. APEX spectra of NGC1333-IRAS 4A smoothed to the 
NGC1333-IRAS 4A and IRAS 4B. Arecibo resolution of 40" (HPBW). 
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Fig. 4. The higher intensity spectrum shows two 290 GHz 
methanol lines from the APEX data towards IRAS 4A that is 
smoothed to 40" resolution and includes contributions from 
IRAS 4B. The weaker spectrum represents the difference be- 
tween that spectrum and another 40" resolution spectrum from 
the same position that does not include data in the direction of 
IRAS 4B. The "contamination" from IRAS 4B is thus estimated 
to be at most 3.5%. 



Fig. 2. 6.7 GHz spectra of NGC1333-IRAS 4A and IRAS 4B 
that have been deconvolved to remove their mutual contami- 
nation. A comparison with Fig. [T] reveals that the spectrum of 
IRAS 4A is almost identical, while IRAS 4B shows a more sig- 
nificant difference. 



shapes of the two lines is striking. From the channel maps of 
the 290.1 GHz line (Fig. |6]l, it can be seen that the emission 
around km s ' and 11 km s ' (LSR) arises from the blue 
and red-shifted lobes respectively of the outflow in IRAS 4A. In 
contrast, the continuum emission from the region at centimeter 
wave lengths is concentrate d over a 0.42" x 0.28" region at the 
core dReipurth et al.ll2002l) . Moreover, the continuum flux den- 



sity at 3.6 cm is 0.43 mJy (a dding the emission fr om IRAS 4A1 
and IRAS 4A2 components : lReipurth et al.ll2002h which is more 
than an order of magnitude weaker than the observed depth of 
absorption (several orders of magnitude weaker when account- 
ing for beam dilution). This indicates that the 6.7 GHz line is be- 
ing absorbed against the cosmic microwave background (CMB) 
over the full spatial extent of thermal methanol emission (includ- 
ing the outflow), as opposed to absorption against the continuum 
from the core. This is similar to the observati on of 12.2 GHz 
absorption against the CMB in dark clouds by IWalmslev et al.l 
(.1988,) . 



4 



J. D. Pandian et al.: Detection of 6.7 GHz methanol absorption towards hot corinos 



^ llllllllllllll 'ftit l ^^ llll|llll| < ll i |llll ^^ llll|l H l|l ^ ll|lll 



o 
o 
o 



31°14'00" 



31°13'00" 



-2.2 "-1.06 " 0.07 



|-u I I I I I I I I I I I I I I I m 

1.21 



+ 



+ 



+ 




+ 




3.48 





ii I 1 1 ^ ^ iii i. |iiii| ku i |tTr i ^^ iiii|iiii | iiii|iii 



9.16 



10.29 



11.43 






ffl I I I I I I I I I I I I I iJ-TTTl ^ I I I I I I I .1 I I I I I I I I I I I l-H I I I I I I I I I I I I I I I ^ I I I I I I I I I I I I I I I I I I I I- 




llll|llll|i m[j lll ^^ l -^ l|llll| U; ll|llll ^ 



12.56 



+ 



+ 



I I I I I I I I I I n I I I I 

2.34 




1 1 1 1 1 1 1 \ \ \\ _ \ Jj \ 1 1 1 1 [ 



13.7 



+ 



+ 



a (J2000) 

Fig. 6. Channel maps of the CH3OH 60 - 5o transition at 290. 1 1 GHz obtained with APEX and smoothed to a spectral resolution 
of 1.1 km s The channe l velocity is written in the top left corner of each panel. The crosses mark the positions of the protostars 
IRAS 4A1, 4A2, and 4B (" Loonev et al]|2000l) . The emission below 4 km s ' and above 1 1 km s ' is associated with the bipolar 
outflow powered by IRAS 4A. The contour levels are at integer multiples of 0.25 K. 



4. Statistical equilibrium calculations 

4. 1 . Anti-inversion in fhe 5 1 - 60 6. 7 GHz transition 

As first realized by lBatrla et al.l (Il987h . methanol maser s come in 
two varieties, termed Class I and II by |Mentenl(ll991 alibi) . Class 
II methanol masers are closely associated with young high-mass 
(proto)stellar objects and arise from the same regions as hy- 
droxyl (OH) masers. Like OH masers. Class II methanol masers 
are most likely pumped by far-infrared radiation in dense (n ~ 
10'^ cm-^), warm (T ~ 150 K) gas Ce saroni & Wal mslevlll99H 
Cragg et al.l 120021: ISobolev & Deguchii il994a.ibi: iSobolev et al.l 
1997). In contrast. Class I methanol masers are frequently found 
in the general vicinity of intermediate and high-mass star forma- 
tion, but often significantly offset (up to a parsec) from promi- 
nent center of acti vity, such as infrare d sources or compact ra- 
dio sources (e.g., iMenten et al.lll986 ). Moreover, their pump- 
ing mechanism can be explaine d from basic propertie s of the 
methanol molecule ( lLeesl ll973t iLees & Haqu'3 11974 . Often, 
Class I methanol masers are accompanied by absorption in the 



Class II maser lines against eith er the CMB (e.g.. lWalmslev et aTl 
11988'; 'W hiteoak & Peng! '19 891) or the continuum of the sourc e 
(WhiteoiietalJll988; Mentei1 ll99l"bl:lPeng & Whiteoakll 19921) . 
The excitation properties of the 12.2 GHz line in absence of 
a strong infrared radi ation field were explained in detail by 
IWalmslev et all (Il988l) . The 6.7 GHz line was predicted to 
be anti-inverted (i.e. over-cooled) under conditions si milar to 
those pr odu cing anti-inversio n in the 12.2 GHz line bv iMentenI 
(1991b') and'Crag g et"aLl ^92). To our knowledge, the 6.7 GHz 
line has been detected in abso rption mostly against the radio con- 
tinuu m of the source (e.g., lMenteniri991bl: llmpelhzzeri et al.l 
I2OO8I) . and the only detections of absorption against the CMB 
are what is repo rted here, and p robably, in the star forming re- 
gion NGC 2264 ( Mentenlll991bl) . 

To investigate the physical conditions leading to absorp- 
tion of the 6.7 GHz line against the CMB, we carried out a 
large velocity gradient (LVG) analysis with spherical geome- 
try (Leuriniet al. 2004, 2007). In contrast to the approach used 
bv lLeurini et al.l (l2007h for modeling high-mass star forming re- 
gions, the only external radiation field used in our calculations is 
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Fig. 5. A comparison of the line profiles of 290.1 GHz emission 
(left axis) and 6.7 GHz absorption (right axis) of methanol in 
NGC1333-IRAS 4A. The location of peak emission/absorption 
is indicated by the dashed line. The similarity of the line profiles 
indicates that the 6.7 GHz line is being absorbed over the full 
spatial extent of thermal methanol emission. 



that of the CMB. Moreover, only a sub-sample of energy levels 
(the first 100 levels) was used for the calculations. 

Figure |2] shows the results of our calculations for the 5i - 
6o transition at 6.7 GHz, for two different CH3OH-A column 
densities (10^"* and 10'^ cm"^), at temperatures of 10, 30 and 
50 K. A linewidth of 1 km s"' was assumed for all cases. For 
comparison, the 2o - 3_i-£' line at 12.2 GHz is also shown. It 
can be seen that both lines have a very similar behavior. They 
are strongly anti-inverted at low densities (n < 10^ cm"^) in ab- 
sence of an infrared radiation field, the primary difference being 
that the 6.7 GHz line is quenched at higher densities compared 
to the 12.2 GHz line. For both lines, the intensity of absorption 
increases with abundance of CH3OH. Moreover, the density at 
which the absorption line intensity is maximum, decreases with 
increasing temperature. Since the only radiation field in our cal- 
culations is that of the CMB, our results are not applicable for 
sources that have an IR radiation field, whose effect among oth- 
ers is to produce strong masers in both lines. For a discussion 
of the excitation of the 6.7 a nd 12.2 GHz in the presenc e of 
IR radiation fields, w e refer to ISobolev & Deguchil ('1994b) and 
ISobolev et al.l (Il997h . or to the more recent work of iCragg et al., 
(I2005h . 

We can verify the predictions of our analysis by simultane- 
ously modeling the CH3OH spectrum at 6.7 GHz and 290 GHz, 
assuming that the lines at the two wavelengths trace the same 
gas. In the next section, we will model the CH3OH emission of 
NGC1333-IRAS 4A. As discussed in section ^ we refrain from 
a similar analysis for IRAS 4B, as the source lies near the edge 
of the area mapped with APEX, preventing us from obtaining a 
spectrum of the source at the spatial resolution of the Arecibo 
data. 



4.2. Methanol around IRAS 4A 

As explained in ^ to model the two datasets, we smoothed the 
290 GHz spectrum to the spatial resolution of the 6.7 GHz ob- 
servations (HPBW of 40"). The similarity of the line profiles at 
6.7 GHz and 290 GHz (Fig. |5]l strongly suggests that the ab- 
sorption at 6.7 GHz arises from the same gas that gives rise to 
290 GHz emission. Hence, to interpret the APEX and Arecibo 
CH3OH data, we applied the method described by Leurini et al.l 
(2004), based on an LVG analysis and simultaneous fitting of 
the entire spectrum with different non-interacting velocity com- 
ponents, each of them being characterized by its source size, ki- 
netic temperature, column density and density. 

The interpretation of our observations is rendered difficult 
by the complexity of the IRAS 4 region and by the poor spatial 
resolution of our data. The 40" Arecibo beam includes the hot 
corino IRAS 4A, and the molecular outflow originating from it 
(Choi 2001, 2005). As discussed in section ^ contamination 
from IRAS 4B at this position is small enough to be neglected. 

The line profiles are characterized by strong blue-shifted and 
weaker red-shifted emission. The same asymmetry between the 
blue- and red-shifted emissi on that can be seen in our CH3OH 
data was observed by Belloche et al.l (12006) in CS. From the 
channel maps at 21.7" resolution of the CH3OH {6k - 5k) 
band (Fig. |6]l, we conclude that the blue-shifted and red-shifted 
emission come from compact sources to the north and south of 
IRAS 4A respectively. By fitting the peak of the integrated in- 
tensity map of the blue-shifted wing of the 60 - 5o-A^ line at 
290.1 GHz (in the range of -2.5 to 4 km s ') with a two di- 
mensional elliptical Gaussian and deconvolving for the APEX 
beam, we estimate a source size of 14" for the blue lobe. A 
similar source size is obtained for the red lobe, using the inte- 
grated emission of the red-shifted wing of the H2CO (4o,4 - 3o,3) 
line, which is stronger than the CH3OH (6^: - 5k) transitions. 
However, the signal to noise ratio in the red wing is lower than 
that of the blue wing, resulting in higher uncertainties. Emission 
around the central velocity is centered close to IRAS 4A, but is 
elongated along the north-south direction and is probably con- 
taminated by the blue and red lobes. This is also seen in the 
CH3OH (6±2 - 5±2)-E lines, which are too weak for a channel 
map analysis, but whose integrated intensity map shows three 
peaks, corresponding to IRAS 4B, and the blue and red lobes of 
the molecular outflow from IRAS 4A. However, the peak corre- 
sponding to the blue lobe of the outflow extends to the core of 
IRAS 4A. 

We modeled the spectrum using three velocity components, 
each having a compact size of 14", corresponding to the blue- 
and red- shifted emission and the hot corino in IRAS 4A. Since 
the red lobe emission is detected in fewer transitions, we mostly 
use this component more to obtain a reasonable fit of the line 
profiles than to derive its physical conditions. Figure|8]shows the 
best fit model spectrum for the 6.7. GHz line and the (6k - 5k) 
band, overlaid on the observed spectra. The best fit parameters 
are: T = 17 K, n = 2 lO*- cm-\ and A?(CH30H) = 4 lO'** cm^^ 
for the core velocity component; T - 27 K, n = 4 10^ cm"^, 
A^(CH30H) - 2 10^^ cnT^ for the blue-shifted component; and 
r = 29 K, « = 1 10'' cm ^ A^(CH30H) = 10'^ cm ^ for the red- 
shifted component. The linewidths used in the model are 1.9, 
7 and 8 km s ' for the core, blue-shifted and red-shifted com- 
ponents, respectively. The total column densities are derived as- 
suming the same abundance for the two symmetry states, A and 
E. 

The 290 GHz observations include the 3o - 2_i E line of 
methanol at 302.369 GHz from the upper sideband of the APEX- 
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Fig. 7. Results of statistical equilibrium calculations for the 2o - 3-i-E (left column) and 5i - Gq-A^ (right column) transitions. The 
solid and dashed lines show the predicted line intensities for CH3OH column densities (for both E and A species) of 10^* cm~^ and 
10'^ cm"^ respectively, and a linewidth of 1 km s 
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2A receiver. However, since the sideband ratio of the receiver as 
a function of frequency is not well known, we did not include 
this line in our calculations. We did compare the line intensity 
predicted by the best fit model with the intensity in our data. The 
model predicts a line intensity Tmb ~ 0.6 K against the observed 
value of 0.2 K. Since the sideband ratio of the receiver at this 
frequency is not known, and since no other solution that fit all 
lines was found, we believe our results to be compatible with the 
data. 

We refrain from a more rigorous analysis of our results due 
to the complexity of the region, coarse resolution of the data, 
and the limited number of transitions available for modeling. 
Our observations are mostly sensitive to the density of the gas 
and the CH3OH column density. However, since the emission is 
optically thin, similar results can be achieved by increasing the 
column density of one of the components and decreasing its cor- 
responding source size and vice versa. Finally, the low tempera- 
ture derived by our analysis is mostly constrained by the weak- 
ness of the 6.7 GHz absorption. Although the fit is not excellent 
(with a reduced of 2), the model predictions for the 6.7 GHz 
absorption are compatible with our observations at 290 GHz. 

4.3. Why do we not observe class II methanol maser 
emission? 

According to ICragg et al.1 (l2005h the 6.7 GHz CH3OH line 
shows strong maser action for dust temperatures. To, above 120 
K. To estimate the radius, r^o, of the region around IRAS 4A 
within which that value is attained, we use the Stefan-Boltzmann 
law, since the dust emission is certainly optically thick at all 
relevant wavelengths. We assume a bolometric l uminosity, L,, 
of 1 1 . 6L0, where we have scaled the estimate of ISandell et al.l 
( Il99lh to the Hipparcos distance of the Per OB2 association, i.e., 
318 pc |de Zeeuw et al. 1999). We calculate rna to be 37 AU. 
Optically thin dus t emission, e.g . with an emissivity index, /3, of 
0.85, derived bv .Belloche et al.1 (12006. see their eq. (1)), would 
yield an even smaller ri2o of 22 AU. Modeling IRAS 4A's cir- 
cumstellar environment as a collapsing envelope. iBelloche et al.l 
(2006) derive relations between the density, n, and radial dis- 
tance, r, for r > 500 AU at diff'erent times. At r = 500 AU, their 
"best fit" model predicts a density of n = 10^ cm"^ and their data 
suggests a density that may even be a factor of 3 greater. While 
extrapolating the density to ri2() (37 AU) may not be straight- 
forward, it is almost certain that the density that close to the 
central object will be vastly in excess of the 10*^ c m~^ at which 
any m aser action would be quenched (see Fig. 3 of ICragg et"aLl 
l2005h . 

This realization may be extended to low-mass protostars/hot 
corinos in general and explains the notorious absence of class 
II methan ol masers around s uch objects (see, a.o, the extensive 
survey of iMinier et al.ll2003 ). In plain words: for low-mass pro- 
tostellar objects, the region in which the radiation field would 
allow maser pumping is so dense that the maser emission is 
quenched, i.e., the 6.7 GHz transition's energy levels are ther- 
malized. 

5. Conclusions 

We have made the first conclusive detection of absorption of 
the 6.7 GHz line against the CMB towards hot corino sources 
NGC1333-IRAS 4A and NGC1333-IRAS 4B. A LVG analysis 
indicates that the 6.7 GHz line is indeed strongly anti-inverted 
for densities lower than 10^ cm"^ almost independent of tem- 
perature. Using observations of submillimeter 6k - 5k lines 



of methanol using the APEX telescope, we model NGC1333- 
IRAS 4A, and are able to reproduce the observed 6.7 GHz ab- 
sorption using cold (T ~ 15 - 30 K), dense (n ~ 10'' cm"-') 
gas. The absence of maser emission at 6.7 GHz in hot corinos, 
and low-mass protostellar objects in general, is presumably due 
to the very high densities in the region where the radiation field 
would provide for maser pumping. 
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Fig. 8. Methanol spectra at 6.7 GHz (top) and 290 GHz (bottom) towards NGCl 333-IRAS 4A. The solid line represents the synthetic 
spectrum for the best fit parameters. The inset in the bottom panel shows the 6_i - 5_i-£ and 6o - Sq-A"^ Unes in greater detail. 



